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Abstract 

We analyse and compare the finite-temperature electroweak phase transition 
properties of classically (non)conformal extensions of the Standard Model. In 
the classically conformal scenarios the breaking of the electroweak symmetry is 
generated radiatively. The models feature new scalars coupled conformally to 
the Higgs sector as well as new fermions. We uncover the parameter space lead¬ 
ing to a first-order phase transition with(out) the Veltman conditions. We also 
discuss dark (matter) aspects of some of the models and compare with existing 
literature when appropriate. We observe that to accommodate both, a first or¬ 
der electroweak phase transition, and a phenomenologically viable dark matter 
candidate requires to go beyond the simplest extensions of the Standard Model. 
Furthermore these extensions must all feature new degrees of freedom that are 
naturally lighter than a TeV and therefore the associated models are testable at 
the upcoming Large Hadron Collider run two experiments. 
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1 Introduction 


The Standard Model (SM) is one of the most successful theories of nature. It is nevertheless 
incomplete and suffers of a numbers of theoretical shortcomings. For example it is unable 
to account for the dark matter (DM) problem and cannot explain baryogengesis. 

It is therefore relevant to investigate minimal extensions that can also address some of 
its flaws. In particular we would like to investigated extensions that can lead to a strong 
first order electroweak phase transition (EWPT). 

We know, in fact, that any model attempting to address the baryon asymmetry in the 
Universe should fulfill the three Sakharov’s conditions [1]: The presence of baryon number 
violating operators; C and CP violations and out of equilibrium conditions. The reasons why 
the SM cannot accommodate a baryon asymmetry is that the associated EWPT is a smooth 
cross over, for the physical Higgs boson mass value [2], Furthermore the SM is not a source 
of sufficiently large CP violation that is needed to generate the observed particlc-antiparticlc 
asymmetry. 

If we insist on having an EW scale induced baryogenesis we must extend the SM in 
such a way that the Sakharov conditions are met. Here we investigate in particular the last 
condition that requires the presence of a sufficiently strong first order EWPT. 

Because of the vast landscape of possible extensions of the SM we concentrate on minimal 
ones. We consider, therefore, extensions featuring new scalars and fermions and explore in 
detail the associated classically conformal limit where the electroweak symmetry breaks 
radiatively at zero temperature. For each model we present first the non-conformal classical 
limit and then move to investigate in detail the classically conformal limit with(out) imposing 
the Veltman conditions (VC). 

We unveil the region in coupling space allowing for a first order EWPT and, when 
relevant, we also discuss the dark (matter) aspects. Not surprisingly we observe that to 
accommodate both, a first order EWPT, and a phenomenologically viable DM candidate 
requires to go beyond the simplest extensions of the SM. Interestingly, however, the exten¬ 
sions that can lead to a first order EWPT must feature new degrees of freedom that are 
naturally lighter than a TeV. This fact renders the models primary candidates to be tested 
at the second run of the Large Hadron Collider (LHC) experiment. Moreover requiring the 
models to be classically conformal and/or to satisfy the one-loop VCs strongly reduces the 
available parameter space of the couplings substantially increasing their productivity. 

To analyse the finite-temperature properties we use the perturbative approach which 
we now briefly summarise. We study the EWPT using the ring improved 1-loop finite 
temperature effective potential given in Eq. 51 of App. A 1 . We then solve for the quantity 
h c /T c . As indication of the possible occurrence of a strong first order phase transition we 
require h c /T c > 1. Here h c = ( h(T c )) ^ 0 is the vacuum expectation value (VEV) of the Higgs 
field, in the EW symmetry broken minimum of the potential, at the critical temperature T c . 
The quantities, h c and T c , are determined from the potential when two degenerate minima 
appear in the potential separated by an energy-density barrier. 

The different degrees of freedom (DoF) that contribute to the effective potential, Eq. 51, 
are given in App. C, D and E for each model separately. In all our calculations we use the 
Landau gauge, and for each model from the SM side we consider corrections from the gauge 
bosons, Z and W, and the top quark. Following common practice (e.g. [4,5]), we neglect the 

Except, in Sec. 4.3.4 where we use the tree level potential augmented only by the leading finite temper¬ 
ature corrections [3]. 
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the Goldstone boson contributions to the effective potential 2 . Similarly, in the conformal 
case we neglect the Higgs boson contribution. 

In Sec. 2 we summarise the EWPT status for the SM. To set the stage, and to acquire 
insight we also consider the EWPT for the SM with imposed VCs and/or its classically 
conformal limit even though these models are ruled out by experiments. 

We then move to Sec. 3 and Sec. 4 where we systematically extend the SM and investigate 
the order of the associated EWPT. When relevant we also discuss the DM properties of the 
models. We finally offer our conclusions in Sec. 5. We provide the bulk of the computations, 
including the expressions for the effective potential at nonzero temperature, for each model 
in several appendices App. A-F. 


2 (Non) conformal Standard Model 

We start our investigation by briefly summarising the EWPT status in the SM. For illustra¬ 
tion and to set the stage for future examples we also investigate the model with imposing 
the Veltman conditions [8] as well as the limit in which the symmetry breaks radiatively via 
a Coleman-Weinberg phenomenon [9]. 

2.1 Standard Model 

The EWPT within the minimal SM has been studied extensively in the past using pertur¬ 
bative methods [10-12] and non-perturbative lattice simulations [2,13-20]. Based on these 
studies, it is well known, that the EWPT in the SM is not first order but a smooth cross over, 
for the physical Higgs boson mass value. More specifically one finds, via lattice simulations, 
that the Higgs mass should be less than about 70 GeV for the phase transition to be first 
order. 

We shall be using the perturbative analysis to provide a first qualitative understanding 
of the phase structure of several extensions of the SM. This analysis becomes robust in the 
limit when the thermal-loop corrections are reliable. This occurs for certain region of the 
parameter space of different theories. In the SM this occurs for light Higgs masses, i.e. less 
than 70 GeV. The break down of the perturbative analysis, for heavier Higgs boson masses in 
the SM, is due to the so-called Linde’s problem [21] which is caused by the gauge bosons that 
are light near the EW symmetric phase. However in models with extra scalars coupled to 
the Higgs via a portal coupling they can provide an important contribution to the transition, 
de facto, improving on the reliability range of the perturbative analysis. We exemplify this 
issue below. 

To elucidate the potential mechanism that can yield a first order nature of the elec- 
troweak phase transition we consider the simple case when the SM effective potential can be 
approximated to be 


VaS(ft, T) « \m 3 s (T)h 2 - ETh 3 + ^A ,„h ‘, (1) 

where m 2 ff (T) = —/r 2 + A (3 g 2 -|- g ' 2 + 4^2 is the thermally corrected Higgs boson 

mass term giving the leading finite temperature term of the potential T 2 h 2 . The bosonic 
degrees of freedom provide the second term to the effective potential, —ETh 3 , where A is a 
coupling constant dependent function. More precisely, if taken into account only the most 

2 To improve the analysis one could use a resummation technic for these states following [6,7]. 
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relevant gauge boson contributions, E = j^(6(g/2 ) 3 + 3 (\/g 2 ~\~ g' 2 / 2) 3 ). This — ETh 3 term 
plays an important role for a possible first order nature of the EWPTThis follows from 
the fact that it is negative and it has a cubic dependence on the background Higgs field. 
Indeed, the interplay of this term with the positive leading term ~ T 2 /r 2 , and the positive 
quartic term jA e ff/z 4 , where A e g(T) « A allows for the formation of two degenerate minima 
in the effective potential. The barrier separating the two vacua hints to a first order nature 
of the transition. And indeed the transition is first order, as long as the Higgs mass is < 70 
GeV [2]. However, the transition is not strong, i.e. h c /T c < 1, for Higgs mass values larger 
than ~ 10 GeV. 

Note that, even though the approximative effective potential, Eq. 1, qualitatively de¬ 
scribes the EWPT at low Higgs masses, it is not a proper analysis. Indeed, if we use Eq. 1 
and the dominant W-boson contributions in the E factor, we get that h c /T c « 2E/X ps 
g rri^y /(tt m 2 H ). However, this approximation is not reliable particularly for heavy Higgs bo¬ 
son masses. This is due to the previously mentioned Linde’s problem, that occurs when the 
gauge boson mass, m ~ g h, is of the order of the magnetic mass, m ~ g 2 T. In this case 
the high temperature expansion parameter, of the form ~ g 2 T/m , becomes of the order of 
unity, and the perturbative approximation breaks down. 

We can circumvent the problem by requiring gh > g 2 T near the transition, which yields 
h/T > g. On the other hand, we had that h c /T c ~ gm-w/m-H■ Thus, the perturbative 
results holds as long as mw/iriH > 1, where as for h c /T c < 1 the perturbative results are less 
reliable. Finally, notice also that, if the daisy diagram resummation 4 , affecting the — ETh 3 
term, is taken into account, it makes the transition slightly weaker. On the other hand if 
higher order corrections are taken into account they tend to make the transition stronger [ 10 ]. 

Because we will be using the perturbative analysis here, other non-perturbative methods 
should also be employed to firmly establish the nature of the phase transition. We have 
checked our 1-loop SM result for h c /T c against the one of Ref. [10] in Figs. 25 and 27, and 
found a very good agreement. 

2.1.1 Veltman conditions 

The scalar boson(s) suffer from the naturality problem expressed, for the SM, in terms of 
the presence of the quadratic divergences for the Higgs boson mass operator. To alleviate 
this problem, one can impose the VCs [ 8 ] for the scalar mass term. In practice this requires 
to set a relation between the model couplings at each order in perturbation theory, so that 
the quadratically divergent piece cancels. 

Using a cutoff regularisation at the one loop level in the SM the quadratically divergent 
piece reads 


5 A 2m 2 H = 


6 A + ^ g 2 + ^ g /2 - 6 y 2 


A 2 

167T 2 


( 2 ) 


3 This same bosonic term provides the 1st order EWPT also in the conformal Higgs portal models studied 
in this work. In these models it is usually the Higgs portal coupling A hs which dominates the E factor and 
thus also limits the so called Linde’s problem. 

4 The daisy diagram resummation procedure is introduced to cure the IR divergences related to the bosonic 
zero modes [10,22]. 
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( 3 ) 


To this order the VC requires S\ 2 m 2 H = 0° yielding 

6 A + | g 2 + ^ g 12 - 6 y 2 = 0 . 

It is well known that this constraint is not satisfied phenomenologically for the SM. Never¬ 
theless for illustration we show its consequences for the EWPT. We keep the gauge boson 
masses, i.e. g and g' coupling values, as measured and use the VC to either predict the wrong 
Higgs boson mass or the wrong top quark mass i.e. the associated coupling value for A or 
for the top Yukawa y t . 

If we keep the correct top mass value, i.e. m t ~ 173 GeV, the VC, Eq. 3 yields mu ~ 314 
GeV. Thus, in this case, since the Higgs boson is very heavy, perturbation theory is not 
reliable, and as previously discussed, lattice results show that the EWPT is a cross over. 

If, on the other hand, we keep the Higgs boson mass as measured rrin ~ 126 GeV the top 
mass gets fixed to m t ~ 96 GeV. Using perturbation theory, we deduce h c /T c ps 0.1, which 
does not suggest the onset of a strong first order EWPT. For such a small ratio, and because 
the Higgs mass is large, i.e. larger than 70 GeV, we resort to lattice field theory that, as in 
the case of the SM, predicts the EWPT to be a cross over. 

To summarize, the VCs do not alter the fate of the EWPT in the SM, i.e. it remains a 
cross over. 


2.2 Conformal Standard Model 


Now we consider EWPT in the classical conformal limit of the SM. This allows us to in¬ 
vestigate how different assumptions on the couplings and structure of the model affect the 
EWPT. 

In this case the EW symmetry breaking is triggered by the CW mechanism [9]. Here the 
tree level potential V(h) = |/r 4 needs to be flat i.e. it needs to fulfil the condition A ~ 0 at 
the electroweak scale. The CW prediction for the one-loop generated Higgs mass is 


m H = 


8tt 2 


tx (3y + 2g 2 g' 2 + y ) + 4A 2 — y\ 
16 


v ew > 


(4) 


where v ew = 246 GeV is the VEV of the Higgs field. For more details see e.g. App. A.l. 

Using the flatness condition, A ~ 0, and assuming the correct top quark mass gives 
m 2 H < 0, which leads to a further instability, that can be avoided by reducing the (Yukawa 
yt) quark mass to smaller unphysical values. By using Eq. 4 with A ~ 0 and floating y t , 
we get that 0 < < 9.8 GeV with top mass varying in the range 78.6 > m t > 0 GeV 

respectively. Actually the lower (upper) bound for the top quark (Higgs boson) mass should 
not be taken strictly because we have neglected all the other fermions. Thus for example by 
fixing m t ~ 10 GeV we have m# ~ 9.8 GeV and h c /T c « 8, which indicates the occurrence 
of a very strong first order phase transition. Increasing the top quark mass value, and thus 
also lowering the Higgs boson mass value, makes the transition even stronger. 

One should notice however, that, for large values of h c /T c , e.g. h c /T c > 5, the energy 
density barrier between the EW symmetric vacuum and one of the broken minima may 
be so large, that the symmetric phase becomes meta-stable and the model becomes not 

5 Here 5^2 indicates that we are only taking into account mass corrections which are proportional to A 2 , 
leaving out the log(A 2 )-terms. Further, it should be understood, here and here after, that all the couplings 
are defined at the same renormalization point (scale) Q , so that e.g. weak coupling g = g(Q ) has the measured 
value at the fixed scale, for example, at the Z boson mass scale Q = niz- 
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viable. Since, however, the potential barrier in the classically conformal scenario is only 
radiatively induced at nonzero temperature, (and absent at one loop at zero temperature) 
the metastable vacuum could decay. In the other models, studied in this work, we restrict 
our attention to the cases h c /T c > 1 but with the ratio which is not too large h c /T c < 4. 
This is typically sufficient to avoid the meta-stability problem. Therefore we do not perform 
in depth meta-stability analysis [24-27]. 

2.2.1 Veltman conditions 

Classical conformality does not address the naturality problem as discussed in [23]. A way 
to alleviate this problem also for the classical conformal scenario is to further impose the 
VC [23], 

By imposing the VC on the conformal SM limit at the one-loop level, Eq. 3 predicts that 
the top (Yukawa) quark mass should be rrit ~ 73 GeV and mu ~ 5 GeV as shown in [23]. 
Because the top is relatively heavy, and the EW symmetry breaking occurs radiatively, the 
EWPT is very strongly first order as indicated by the h c /T c ps 16 value. We remind the 
reader of a potential meta-stability issue [24-27]. 

To summarize, the SM does not support a 1st order EWPT with or without VC. The 
situation changes dramatically in the case of the classically conformal SM where a strong 
first order phase transition occurs with and without VC albeit the model is not phenomeno¬ 
logically viable. 

Thus, we would like to investigate and compare (classically conformal) extensions of the 
SM with(out) VC that can support a 1st order EWPT. 

A reasonable strategy is to search for models that enhance the negative cubic term 
—ETh 3 . As we shall see this mechanism is particularly efficient for the conformal models 
since they have a flat direction for the effective potential at the tree-level. The enhancement 
is then achieved by adding new bosonic degrees of freedom to the model. This guides us 
towards theories containing new scalar degrees of freedom to be added to the SM. 


3 (Non) conformal Extension with A Real Singlet 
Scalar 

We begin the investigation of extensions of the SM that can support a first order EWPT 
with Higgs portal models featuring a real singlet scalar. This is a heavily studied model and 
it is reported here mostly for completeness. We will refer to the existing literature when 
needed while concentrating particularly on the classical conformal limit with(out) VC. 

3.1 SM + Real Singlet Scalar 

The most general renormalizable tree level scalar potential for this model reads: 

Vo = -^ 2 H ] H + \ h {H ] H ) 2 + ‘ n ^S 2 + \ hs H ] HS 2 + ^S A + 

+ ^S 3 + f 4S + f i m SH^H , (5) 

where S is the new real singlet scalar. The last three terms are absent if one requires the 
potential to be Z 2 symmetric. The hnite-temperature phase structure can be studied by 
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replacing the fields in the potential by their background values h and s. The associated tree 
level potential is obtained by substituting H —» h/\J 2 and S —>■ s in Eq. 5. Depending on 
the values of the model parameters the phase structure and the vacuum structure can be 
very rich. We shall not make a full EWPT analysis for this general model here, as it can be 
found from the literature [3,28-32] (see also [33]). A thorough investigation of the potential 
structure, yielding a strong first order phase transition stemming already from the tree level 
potential, once the leading finite temperature terms are taken into account, was recently 
performed in [3]. We will consider this kind of scenario later in Sec. 4.3.4. 

The EW symmetry breaking and first order phase transition can also be obtained via 
loop corrections. This is argued in Sec. 3.2 where we will discuss the EWPT for the classical 
conformal case with(out) VC. 

If the Z 2 symmetry is imposed the model can also feature a DM candidate because the 
singlet becomes stable. The DM phenomenology and EWPT of this type scenario has been 
studied recently e.g. in [34,35]. The authors concluded that in order to have a sufficiently 
strong first order EWPT capable to produce baryogenesis, the singlet can only be a sub¬ 
dominant DM component. 

3.2 Conformal SM + Real Singlet Scalar 

In the classically conformal case the tree level scalar potential reads: 

Vo = A h (H'H) 2 + \ hs H'HS 2 + jS 4 , (6) 

which automatically respects the ^-symmetry, as no terms with dimensional couplings, 
including the terms with odd number of S fields, are allowed by the classical conformality 
requirement 6 . Stability considerations constrain the model parameters to be 

Xh > 0 and A s > 0 , and if Ah s < 0 , then Aj[ s < AhA s . (7) 

In this case the EW symmetry breaking follows from the CW mechanism. Because the CW 
mechanism, to be operative, requires the presence of flat directions at the electroweak scale 
we can consider different scenarios. Here we shall consider two cases. The case (1) in which 
the coupling A h s is positive and the flat direction is in the /i-direction. In this case S does 
not acquire a non-zero VEV, i.e. (s) = 0 at zero temperature (and at finite temperature). 
In case (2), A hs is negative, and the flat direction is in some non-trivial direction in the 
(h,s)-plane. Here the field S acquires a non-zero VEV, i.e. (s) 7 ^ 0 at zero temperature (and 
at finite temperature). 

3.2.1 Case (1): (s) = 0 

For A hs positive and with the scalar quartic couplings satisfying the relation A s > A^ = 0, 
the S field does not acquire a VEV and the flat direction can be taken to be in the Higgs 
direction. 

6 Of course classical conformality per se does not always imply the presence of a ^-symmetry. This is 
apparent from the case studied in Sec. 3.4, where a new singlet Weyl fermion, with a classically conformal tree 
level interaction term y x Sxx + h.c., is added to the model given in Eq. 6. Therefore unless a Z 2 -symmetry is 
strictly enforced, the scalar S does not need to be stable, even though the tree level scalar potential possesses 
an accidental ^-symmetry. Furthermore if S acquires a VEV one should expect the presence of new sectors 
breaking explicitly the .^-symmetry to avoid the domain wall problem [36,37]. 



2.0 


1.8 

1.6 

hc/T c i.4 
1.2 

1.0 

0 1 2 3 4 5 6 

A s 

Figure 1: Plot of h c /T c as a function of \ s for the real singlet Higgs portal model of case (1). 
Here the Higgs boson has mass 126 GeV and the singlet S has mass 541 GeV. 


The phase transition for this type of model has been investigated previously in [5,38]. 
There, however, the singlet quartic self coupling effect was not included but on the other 
hand the model included several real singlet scalars with identical portal couplings A hs- We 
were able to reproduce their results and find that in the case of 12 real singlet scalars, 
the EWPT can be strongly Erst order with h c /T c « 2.7, provided the portal coupling is 
A hs ~ 1.4 7 . 

We used the full one-loop effective potential, Eq. 51, with mass parameters, as defined 
in App. D. The flat direction condition is A h(Q) ~ 0, with Q ~ vew , and the GW analysis 
is similar to the one for the SM but able to provide the correct Higgs boson mass ~ 126 
GeV. In the case of only one singlet scalar, and neglecting A. s , we find that h c /T c ps 1.76 for 
A hs ~ 4.84 which also allow for the correct Higgs boson mass. The singlet tree level mass, 
follows from the relation m 2 = A h s v\ w and yields m s ~ 541 GeV. 

We show now how the inclusion of A s affects the above results. The one-loop Higgs mass 
and the tree level singlet masses are unaffected by the addition of the S quartic interaction 
term. A s affects the EWPT because it appears in the singlet scalar thermal mass. Increasing 
A s decreases the ratio h c /T c . This follows from the fact that the temperature dependent 
piece, in the daisy improved boson-induced cubic term in the effective potential increases 
with A s . This damps the held dependent piece responsible for the 1st order nature of the 
transition. In turn, one observes a decrease of the ratio h c /T C) as it can be seen from Fig. 1. 

Finally, using the DM results from [34], we conclude that in the conformal case the singlet 
can provide at most a subdominant DM component with density < 1% of the observed DM 
density. 

For the Higgs portal models we allow relatively large couplings, i.e. up to < 47T. Of 
course, if some of the parameters approach 47 t it is clear that the one-loop approximation 
is not reliable. In this case one should take into account higher loop effects, or better use 
directly non-perturbative methods. However, recently in [39], for complex singlet Higgs 
portal models it has been shown, by an explicit 2-loop computation, that identifying a 1st 
order EWPT via the 1-loop effective potential provides reliable results up to quite large 
values of the couplings of the order of A hs ~ 8. 

7 The notation for the portal coupling in [5,38] is different from ours, the two are related via A hs = £ 2 - 
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3.2.2 Veltman conditions 

Now we impose the VC for this model. In this case we have two conditions, one for each 
scalar held. The singlet scalar VC is defined similarly as the Higgs one, i.e. 8{^m 2 s = 0. The 
conditions read for S and H respectively [23] 

3A S + 4A/j S = 0 , (8) 

6A h + ~^g 2 + ~^g' 2 ~ 6 y 2 + A h s — 0. (9) 

We observe immediately from Eq. 8 that A hs should be negative since A s > 0 to ensure 
stability of the potential. However, the EWPT study and CW analysis, in Sec. 3.2.1 above, 
indicated, that A hs should be positive, to be able to support a strong first order EWPT and 
the correct Higgs boson mass value. Thus the VCs are naturally applicable to case (2) rather 
than (1) for which A^ s < 0 from the start. 

3.2.3 Case (2): (s) ^ 0 

The hat direction is not in the Higgs direction but it is chosen to be along some non-trivial 
direction in the (h,s)-plane. For the CW analysis we follow [23] and review the hat direction 
conditions in App. A.l. 

It is convenient to write the potential, Eq. 6, using polar coordinates, i.e. H = 
(0, r cos cu) T /\/2 and S = r sin u which now reads: 

^4 

Vo = — (Ah cos 4 o; + A s sin 4 a; + 2\ hs sin 2 u> cos 2 u>) . (10) 

The minimum of the potential, along some ray r, is found by studying the hrst and second 
derivatives of the potential, leading to the conditions 

0 < X h < min{A s , X hs } : (u>) = 0 , (11) 

0 < A s < minjAft, X hs } : (u) = ^ , (12) 

~VX h X s < x hs < minjA/j, A s } : tan 2 (o;) = ^^ . (13) 

As Afis 

We chose to use polar coordinates, because is it very useful for the CW and EWPT analysis, 
as we will now discuss. Indeed, once the model parameters X s , Xh and A hs are such, that 
the hat direction condition is fulfilled, then the hxed angle ( oj ) defines simultaneously the 
minimum direction and the hat direction of the tree level potential, while (r) the location of 
the ground state along the hat direction. Thus once this direction is hxed to (c u), the CW 
and EWPT analyses can be performed directly along the ray r, see App. A.l for a detailed 
discussion of the method. 

Now, the hrst case, (u>) = 0, given in Eq. 11 above, corresponds to the case (1) studied 
in Sec. 3.2.1 above. On the other hand for ( oj ) = |, given in Eq. 12, the EW symmetry does 
not break and we are left with the third possibility, i.e. Eq. 13. The hat direction condition 
in this case reads 


XhXs — X hs — 0 . 


(14) 
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The CW mechanism along the flat direction, Eqs. 14 and 13, leads to spontaneous symmetry 
breaking with (r) cos(u;) = Vew■ The flatness and VCs are simultaneously satisfied at 
Q ~ Vew■ By expanding our fields around the ground state via rcos(o;) = v ew + h 0 and 
rsin(a;) = vew tan (a;) + s 0 , we find the mass eigenstate fields 


(ft = ho cos (u>) + so sin (a;) , $ = so cos((n) — ho sin(a;), 


(15) 


with tree level masses = 0 and = 2(A h — Xh s )v 2 EW . Here 0 is the physical Higgs 
field and 4> the heavier state. The CW induced one-loop mass for the Higgs field is 


ml = 


cos 2 (a;) 
87T 2 


i(6 9 4 + 3( 9 2 + g' 2 ) 2 ) + 4(A„ - \ H ,) 2 - ^y? 


- EW 


(16) 


We can now constrain two out of the three model parameters Xh, X s and A h S , using the flat 
direction condition, Eq. 14, and the Higgs boson mass relation, Eq. 16. We choose to first fix 
A hs using the flat direction condition, Eq. 14. By taking into account the stability conditions 
A h, s > 0 and the minimum condition, Eq. 13, this gives A hs = ~%/XhK- Then, using the 
previous relation, as well as the measured y t ,g and g' coupling values, in the Eq. 16, we 
can search for real valued solutions for A h, as a function of A s , providing the physical Higgs 
boson mass value. Given that A h s and A h are now functions of A s , we then performed a scan 
over the last remaining free parameter A s , in our EWPT study within the range 


0.01 < A s < 47T, 


(17) 


with a flat prior distribution with 1000 points. In the phase transition analysis we used the 
1-loop effective potential, Eq. 51, with the relevant mass terms given in App. D and with 
the CW conditions discussed above and given also in App. A.l. 

In Fig. 2 we show how different model parameters relates to X s , whereas the distributions 
of the model parameters are shown in Fig. 3. Same data is used in both figures and all shown 
points provide h c /T c > 1. In the lower(middle) left panel of Fig. 2 (Fig. 3) we show how 
the m 0j $, i.e. the mass of the new scalar, relates to A s (is distributed). The parameter 
cos(cu), in the lower right panels of Figs. 2 and 3, determines how much the Higgs couplings 
change, when compared to the SM values. The SM couplings are reproduced for cos (a;) = 1. 
The purple points/histograms show the distributions once the constraint cos(cj) > 0.85 is 
taken into account. The value cos(a;) = 0.85 is the 2cr lower bound from the best fit value 
cos(cu) = 0.95 to the experimental results [40]. All the points are for the correct Higgs boson 
mass ~ 126 GeV. From the top most middle panel of Fig. 3 we note that for the model to be 
consistent with cos(cj) > 0.85 the singlet quartic self coupling X s should assume rather large 
values. For the smallest values, A s ~ 5, the one-loop perturbative approach might still be 
a reasonable approximation. Notice also, that $ cannot be DM here because it is unstable. 
Indeed, $ is heavy enough (see lower left panel of Fig. 2) to decay into two Higgs fields. 

To summarise, the model can provide a strong 1st order EWPT and correct Higgs boson 
mass. However, to avoid the experimental constraints set for cos(cj), the singlet quartic self 
coupling tends to have rather large values X s > 5. The DM problem cannot be addressed 
with this model. 
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Figure 2: Real singlet Higgs portal model parameters as a function of A s for case (2) that 
yield a first order EW phase transition. The purple (blue) points with(out) the constraint 
cos(cu) > 0.85. For all the model points the Higgs boson mass assumes the correct value. 
The model parameters have not been constrained using the VCs. 



3.2.4 Veltman conditions 

The Veltman conditions are Eqs. 8 and 9. Once the VC condition, A h s < 0, from Eq. 8 is 
satisfied the flat direction condition, Eq. 14 reads 

V A/iAs + A hs — 0 • (18) 

This condition agrees with the flatness condition in Eq. 18. Using therefore the VCs, Eqs. 8 
and 9, along with the flat direction condition, Eq. 18, the model parameters get constrained 
via 


3 16 4 Q 

A*, = --A s , A s = -A* , cos 2 (a;) = - , X h — -(8 y 2 - 3 g 2 - g' 2 ). (19) 

Implementing these conditions on the Higgs boson mass relation Eq. 16, and using the 
measured values of the top Yukawa and the SM gauge couplings, we get 

~ 95 GeV and m$ « 541 GeV. (20) 

Now all the model parameters are fixed and using these values we End that h c /T c « 2.2 
indicating the occurrence of a strong Erst order EWPT. Because the mass of the lightest 
scalar is lower than the one used in the Fig. 3 this model is not present there. Although the 
resulting Higgs mass does not agree with experiments to this order in perturbation theory, 
it would be interesting to investigate higher orders to see if they might alleviate the tension. 

To summarize, once the one-loop VCs are imposed, and holding fixed the top mass to the 
experimental value, we predict a Higgs boson mass of ~ 95 GeV and strongly first order 
EWPT as indicated by the h c /T c « 2.2 ratio. The SM particle couplings to the physical 
Higgs boson are suppressed by a factor cos(o;) compared to the SM. More precisely the model 
predicts cos(o;) ~ 0.76, which is below the experimental 2cr-limit cos((u) = 0.85 [40]. 
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Figure 3: Real singlet Higgs portal model parameters distributions for case (2) that yield 
a first order EW phase transition. The purple (blue) histograms show the distributions 
with(out) the constraint cos(cu) > 0.85. For all the model points the Higgs boson mass 
assumes the correct value. The model parameters have not been constrained using the VCs. 


We learnt that the classically conformal singlet Higgs portal model can support a strong 
1st order EWPT and provides the correct EW spectrum. This can be viewed as an improve¬ 
ment on the SM. Further imposing the VC reduces the Higgs mass but strengthen the 1st 
order EWPT. Finally, the model can provide at most a subdominant DM component. 

We move now to the next natural step, i.e. adding the simplest new fermionic sector to 
the theory. 

3.3 SM + Real Singlet Scalar + Weyl Fermion(s) 

Now we shall consider a model, where the simplest singlet scalar Higgs portal model is extend 
with a singlet Weyl fermion y. The general non-conformal scalar potential, including the 
new Weyl fermion y interaction term reads 

Vo = -fJ 2 H^H + X h (H^H ) 2 + m 2 s S 2 + X hs H^HS 2 + ^S 4 (21) 

+ fi 3 S 3 + /ilS + + (y x S + M x ){XX + XX) ■ 

The EWPT for this kind of models where, however, the new fermion y is a Dirac particle 
has been studied recently in [40-42], There the particle y was associated with dark matter, 
and it is was found, that the EWPT can be strongly first order and that the model can yield 
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the correct DM density 8 . Here too x can be associated with dark matter [43] and since the 
new fermion sector has a milder effect than the new boson one on the phase transition we 
expect a strong first order phase transition [40-42], 

We will therefore move to investigate the classically conformal limit of the model. 


3.4 PNC model: Conformal SM + Real Singlet Scalar + Weyl 
Fermion(s) 

Once we take the classical conformal limit of the Eq. 21, the tree level potential reduces to 

Ho = A h (H f H) 2 + A hs H'HS 2 + ^S 4 + y x S(xx + XX) ■ (22) 

This model resulted to be the most intriguing among the class of the perturbative natural 
conformality models (PNC) investigated in [23]. The PNC models attempt to provide the 
correct Higgs boson mass via the CW mechanism and to alleviate simultaneously the nat- 
urality problem by satisfying VC. Obviously these models are much more constrained than 
either of the two conditions satisfied independently. The above model was, indeed, found to 
be an ideal candidate PNC model as we shall briefly review below. It is therefore relevant 
to see if the model is also capable of leading to a strong 1st order EWPT. 

We also observe that the new classically conformal interactions of S with the Weyl fermion 
directly destroys the accidental Z 2 symmetry of the scalar part of the tree potential. 

We shall consider two cases: In case (1) only H acquires a non-vanishing VEV. Where 
as in case (2) S can also get a non-vanishing VEV. 


3.4.1 Case (1): (s) = 0 

Here the EWPT investigation is very similar to the one in Sec. 3.2.1. The difference resides 
in the fact that the singlet dark fermion affects the singlet scalar thermal mass. The relevant 
mass terms for the EWPT analysis are given in App. D.3.1 and D.l. 

As further constrain we assume the physical Higgs boson mass which allows, as we shall 
momentarily see, to constrain Xh s . We also assumed that the flat direction occurs along the 
/i-direction and thus we have A ~ 0 at the EW scale. The CW analysis then yields the 
one-loop mass for the Higgs boson 




f(2g 4 + (. g 2 + g' 2 ) 2 ) + \\l 


vt 


J EW 


(23) 


Knowledge of the physical Higgs boson mass value as well as the SM couplings y t ,g and g', 
yields A hs ~ 4.84. We are left with two free model parameters that we allow to range in 


0.01 < A s < 4.5 and 0.01 < y x < 6.5 , (24) 

and search for the parameter space where h c /T c > 1. Our results are reported in Fig. 4. 
Similarly to Sec. 3.2.1 we observe that the extra contributions to the singlet scalar thermal 
mass decreases the ratio h c /T c and as expected the largest values for h c /T c is obtained when 
X s and y x are set to zero. From the Fig. 4 we see that the model can provide h c /T c > 1 
values for a large part of the model parameter space, particularly favouring small A. s and y x 
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Figure 4: We show the h c /T c contour plot as function of the dark Yukawa coupling y x and 
the scalar S quartic self coupling A s . The portal coupling A hs is fixed to the value ~ 4.84, 
which provides the correct Higgs boson mass and satisfies the VC for the Higgs boson mass 
term. On the upward turning blue curve also the VC for the new scalar boson S mass is 
satisfied. 


values. Further, using X hs « 4.84, the model predicts that the tree level singlet scalar mass 
is 541 GeV, as it is the case in Sec. 3.2.1. 

3.4.2 Veltman conditions 

The VCs read [23] 


3A S + 4A hs — 8 y x — 0 , 


(25) 


for the singlet, and 


6A h + ^g 2 + ^g' 2 - 6 y 2 + X hs = 0 , (26) 

for the Higgs field. Compared to the case studied in Sec. 3.2 in Eq. 8, the VC for the 
singlet has been now modified because of the presence of the dark Yukawa coupling y x . 
This allows A hs to be positive here. Further, note that the VC for the Higgs boson, Eq. 26, 
gives A hs ~ 4.84, which is very surprisingly exactly the value needed to yield the physical 
Higgs mass in Eq. 23. Using this value of A hs ~ 4.84 in Eq. 25 the VC returns the relation 
X s ~ | y 2 x — 6.45, shown as the upward turning blue curve in Fig. 4. Thus, on this curve the 
Veltman condition for the scalar boson S is satisfied. To collect, the model predictions are 

O 

m l j, « 126 GeV, « 541 GeV, Ah s ~ 4.84 and A s « -y 2 — 6.45, (27) 

3 * 

8 Notice, the hard mass term M x x\ extends the model parameter space, and plays a significant role in 
the DM studies. 
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that also give h c /T c > 1 once y x <1.7. 

To summarize once the VCs are imposed, the model becomes very constrained and thus 
highly predictive. Remarkably it yields the correct mass for the Higgs boson and predicts 
the mass for the singlet scalar while satisfying the VCs [23]. In addition, the model supports 
a first order EWPT. Although the value of A hs ~ 4.84 is not very small it is still not un¬ 
reasonably large for a perturbative study to hold [39]. Nevertheless higher order corrections 
would be interesting to investigate. Further, the VC for the singlet can be satisfied for very 
small couplings with the limiting case of \ s ~ 0 once y x ~ 1.55 providing also the largest 
ratio h c /T c « 1.4. 


3.4.3 Dark aspects 

Here we briefly discuss aspects related to the new Weyl fermion y. As shown in the previous 
section, once we take the classical conformality and the VC to guide the model building, 
the dark Yukawa coupling y x gets connected to A s in Eq. 27. If we furthermore ask for the 
largest possible h c /T c ratio, the model predicts y x ~ 1.55. Because, however, S does not 
acquire a VEV this Weyl fermion y remains massless at zero temperature. Unless we further 
extend the model y would behave as a sterile light neutrino. If we allow the introduction 
of a small conformal breaking by intruding a hard mass term (or yet another scalar held 
acquiring a VEV) for y then y could be considered as a DM candidate [43]. We will not 
speculate further on this point and move to the next natural case in which also S can receive 
a non-vanishing VEV via the CW mechanism. In this case y will acquire a tree level mass 
term while abiding the classical conformality requirement. 


3.4.4 Case (2): (s) ^ 0 


We will now consider the general case in which also the singlet S acquires a non-vanishing 
VEV. For this case we follow the steps in Sec. 3.2.3. To have a non-trivial hat direction and 
EW symmetry breaking via CW the model parameters should satisfy the relations in Eq. 13 
and Eq. 14. For potential stability reasons the scalar couplings should further satisfy the 
relation A/^ > 0. 

The mostly singlet state has a tree level mass = 2(A h — A hs)v 2 E w as in Sec. 3.2.3. 
However, because the S held receives a VEV, and H and S mix, the one-loop Higgs mass 
gets modified to, 


ml = 


cos 2 (to) 

8lT 2 


^(% 4 + 3 (g 2 + g ,2 ) 2 ) + 4(A h ~ X hs ) 2 ~ ^ yf ~ S2y 4 tan 4 (uj) 


'EW 


(28) 


This means that also y contributes to the Higgs boson mass, even though the effect is 
suppressed by the mixing. Using the hatness condition, Eq. 14, and the Higgs mass relation, 
Eq. 28, we may relate and constrain the model parameters A., A s , XI and y x as in Sec. 3.2.3. 
However, compared to Sec. 3.2.3, here we have one more parameter y x , and thus we made 
the scan over two parameters within the ranges 


0.01 < y x < 5 and 0.01 < A s < 47 t. (29) 

We used hat prior distributions with 2000 points in the scan. In Fig. 5 we show h c /T c > 1 
constant contours in (y x , X s )~ plane. We observe that it is possible to achieve the phenomeno¬ 
logically desirable correct value of the Higgs mass for a very wide range of the parameter 
space along with a 1st order EWPT. 
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Figure 5: We provide the h c /T c contour plot in the (y x , A s )-plane for the PNC model of case 
(2). The experimentally observed value of the Higgs boson mass ~ 126 GeV is enforced. 
The black dot shows the point where the model satisfies also the VCs and it leads to the 
values cos(u;) ~ 0.79, m x = 191 GeV and = 604 GeV. 


3.4.5 Veltman conditions 

The VCs for this model are Eqs. 25 and 26. Using these with the flat direction condition, 
Eq. 14, and the Higgs mass relation, Eq. 23 and using the measured values for the couplings 
g,g' and yt, we can solve for all the four free parameters of the model. Keeping only 
the solutions yielding positive real valued physical quantities, the model parameters are 
constrained to be 


Xh ~ 1.12, X hs « —1.90, A s « 3.20, y x « 0.50, (30) 

yielding the following physical predictions 

TTifj) ~ 126 GeV, m$« 604GeV, m x « 191 GeV, cos(cu) « 0.79. (31) 

Further, using the parameter values Eq. 30 for the EWPT study we find h c /T c ss 1.88 
indicating a sufficiently strong first order phase transition. 

Let’s pause and reflect on the results. Notably, all the couplings are fixed and relatively 
small, thus the perturbative analysis for the phase transition is robust. The model yields 
the correct Higgs boson mass values and predicts one extra singlet scalar and Weyl fermion 
with masses ~ 604 GeV and m x ~ 191 GeV respectively. The Higgs coupling factor 
cos(cu) ~ 0.79 is slightly below the experimental 2cr-limit and higher loop corrections to the 
Higgs mass relation and VCs could further alleviate the tension. Since the H and S fields 
mix we have that x couples to the Higgs via the portal coupling. Thus, if \ is assumed to be 
a DM candidate, because y x is not that small and furthermore m x is naturally predicted to 
be in the classic WIMP mass region, it could be probed via DM direct detection experiments. 
On the other hand \ could equally be associated with a massive right handed sterile neutrino, 
which could then mass mix with the SM neutrinos. Of course, as we have shown the VCs 
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Figure 6 : Distributions of the PNC model parameters in case (2) producing first order 
EW phase transition. The purple (blue) histograms show the distributions with(out) the 
constraints cos(cu) > 0.85 and BR 0 _j. xx < 0.17. The red points, providing a subdominant 
DM component, survive also the LUX [46] constraint. 


can be relaxed an a much more general parameter space is allowed consistent with a 1 st 
order EWPT and phenomenological constraints. 

It is therefore relevant to discuss next the intriguing dark aspects of this model. 

3.4.6 Dark aspects and constraints 

We discuss now shortly about the constraints and aspects related to two different x scenarios. 
The one in which x is the DM candidate and the other in which x is a right handed neutrino. 
We first consider the most evident constraints for the former case, which we call scenario A. 
These constraints can be partially used also for the latter case, called scenario B. 

Scenario A. Here we star by assuming, that x is a DM particle stabilised by a new 
discrete symmetry e.g. Z- 2 . We compute its annihilation cross sections into SM fermions, 
gauge and Higgs bosons as well as into the new scalar 4>. We evaluate the relic density 
using the classic approximate method to check if it can account for the full cosmic DM 
relic density Ddm- The details can be found in the App. F. However, we consider a wider 
parameter space where x can also be a subdominant component quantified via f re i = • 

Here Q x h 2 is the x density parameter (times the dimensionless Hubble parameter squared) 
and Ddm h 2 = 0.1188 is the DM one consistent with the latest observations [44]. Using the 
DM direct detection constraints and the Higgs invisible decay branching ratio bound we can 
now limit the parameter space of the model. 

We first consider the more general model scenario, where the VCs are not assumed. The 
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spin independent y-nucleon cross section for the model is 9 


2 


osi 


4/V 


y 2 x sin 2 (o;) cos 2 (a;) 




( 32 ) 


where p is the y-nucleon reduced mass and / is the effective Higgs-nucleon coupling, for which 
we used the value / = 0.345 [34], To ease the comparison with the DM direct detection 
constraints also for the subdominant case we define the following effective scattering cross 
section (see e.g. [34,35,45]) 

o£ = OSl/rel < C P , (33) 


where (Jgj P is the experimental cross section bound. Naturally, the usual constraint for the 
cross section is obtained once / re i = 1. To constrain the model we use the recent LUX 
results [46] for (Xgj P . 

Let us now explain why it is convenient to use Ogf. Under the expectation that DM 
clusters similarly in all cosmic structures irrespectively on its sub- or dominant nature, the 
factor / re i can be used to properly scale the local y density via the relation p l ° c = freiPvM, 
where p^M ~ 0.3GeV/cm 3 . This affects the interpretation of the DM direct detection 
bounds that depends on the assumed local DM density. Indeed, as the local density for the 
subdominant DM is smaller than the usual p^M> the scattering cross section ogi in Eq. 32, 
can be larger than (Tgj P obtained using / re i = 1. Thus, better quantity to constrain is the 
effective scattering cross section cr|f, defined in Eq. 33, which takes properly into account 
the local DM density scaling for each type DM component. 

In the Fig. 7 we show the scatter plot of the model cross section crgf and the LUX 90% CL 
limit [46] . Colorcoding indicates the relative dark matter density with red giving the full DM 
density where as blue only a percent level fraction. Grey points yield too much dark matter 
i.e. / re i > 1 and thus are excluded. Triangles pass the other experimental constraints, like the 
invisible Higgs decay BR^ XX < 0.17 constraint discussed below, and also the experimental 
2cr-limit cos(o;) > 0.85 mentioned previously. Dots do not pass these extra constraints. 

In the case y is light, such that 2 m x < m^, it can open up an invisible Higgs decay 
channel. In this case the invisible Higgs decay width to 2y is 


T 


'-►xx — 


m# 

47T 


y 2 x sm ! (u) 



(34) 


where is the physical Higgs mass and m x = 2y x t&n(uj)vEw- The invisible Higgs decay 
branching fraction in this case reads 


BID 


0->xx 


►XX 


SM 


+ T. 


<M-xx 


(35) 


where Tsm is the Higgs decay width in the SM. This branching fraction is constrained by the 
LHC data, and for the limit we used BR^ XX < 0.17, where the 95% CL limit BR inv < 0.17 
was obtained from [47]. 

9 The minus sign in the last squared term in Eg. 32 follows directly from the real orthogonal mixing matrix 
which diagonalizes the (h,s) mass matrix. In our case one of the mass eigenvalues is positive and the other 
is zero, i.e. the Higgs-like state has zero mass at tree level. 
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Figure 7: The model effective spin independent y-nucleon cross section <jg® as a function of 
DM mass m x . Colorcoding depicts the relative DM density f re \ for each point. Grey points 
have / re i > 1 and are thus excluded. Large triangles pass the invisible Higgs decay branching 
fraction limit and the experimental 2<7-limit cos(cu) > 0.85, whereas small dots do not. The 
black dot indicates the parameters space point satisfying the VCs. The red solid curve is the 
LUX 90% CL limit [46] above which the model points are excluded. The model data match 
the ones used in Fig. 6 and furthermore for all the shown points h c /T c > 1. 


Let us now check the scenario where the VCs for the model are satisfied. In this case all 
the model parameters are fixed and given in Eqs. 30 and 31. Using these values for the relic 
density analysis we get that fl x h 2 ~ 0.164 providing / rel ~ 1.377, which is surprisingly close 
to the full cosmic DM density. However, once we calculate the scattering cross section, using 
the values from Eqs. 30 and 31 in Eq. 32, we get that ogi ~ 1.56 x 10~ 43 and comparing 
with LUX limit [46], ^(wdm ~ 190 GeV) ~ 2.2 x 10“ 45 cm 2 , we notice that the model 
point, indicated as a black dot in Fig. 7, is ruled out. In fact, this large <7si cross section is 
excluded also by several other direct DM detection experiments. 

To summarize, the model can yield a phenomenologically viable subdominant DM com¬ 
ponent, with < few % of the total observed DM density without the VCs. 

Scenario B. If we associate y, not with DM but simply with a massive right handed 
neutrino we have the same constraints except the DM direct detection limit Osi limit. Thus 
all the purple points in Figs. 6 and 7 are phenomenologically acceptable in this scenario. 
However, once VCs are imposed the model is in tension with the experiments because 
cos(cu) ~ 0.79. 

To summarise, the model allows for the observed value of the Higgs boson mass and first 
order EWPT. It also predicts a new massive scalar and massive Majorana particle. The 
VCs are still very constraining and although the correct Higgs boson mass is still viable 
together with strong EWPT, the Higgs couplings to the SM fields are suppressed by a factor 
cos(cu) ~ 0.79 and deviating by more than two sigma from the LHC lower bound result 
cos (uj) > 0.85. 

We move now to complex scalar extensions of the SM. 


4 (Non) conformal Extension with A Complex Scalar 

In this section we consider the EWPT in minimal complex scalar extensions of the SM. We 
start by briefly discussing the EW singlet complex scalar Higgs portal models, and end with 
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the ”SEyy” model, where the new scalar is not a SM singlet. 


4.1 SM + Complex Singlet Scalar 

The tree level scalar potential of the complex scalar Higgs portal model reads 

V 0 = ~^ 2 H ] H + A h (H^H) 2 + m 2 s S^S + A hs H^HS^S + A s (SS^) 2 . (36) 

Phase transitions in this type of models have been investigated extensively in the past e.g. [39, 
48, 49] 10 . For this reason we shall not make a full EWPT analysis here but concentrate 
instead on the classically conformal limit, and on extensions, that have not been considered 
in the literature. 

4.1.1 Conformal SM + Complex Singlet Scalar 

The classically conformal limit of the potential in Eq. 36 is 

Ho = A h (H'H) 2 + A hsH'H&S + A s {SS*) 2 . (37) 

If only the Higgs field direction is allowed to undergo a radiative CW symmetry breaking 
scenario the EWPT analysis is similar to the one in Sec. 3.2.1. Now instead of one real scalar 
we have two, and adding more real scalars (up to 12) leads to the case studied in [5,38]. 
These results are verified and further improved, for the complex S scalars, via a 2-loop RG 
improved EWPT analysis in [39]. There it is demonstrated, that a strong 1st order phase 
transition can be obtained reliably from perturbative analysis even with large portal coupling 
values. Since this model has been extensively studied elsewhere, we will not repeat here the 
EWPT study. Finally, adopting the results of [34], we note that the model can provide a 
subdominant DM component. 

4.2 SM + Complex Singlet Scalar + Fermions 

If new singlet Weyl fermions Xi are added to the model, one arrives to the Majoron-like 
model [53]. The tree-potential is 

Ho = -fi 2 H^H + X h {H^H) 2 + m 2 s S^S + X hs H^HS^S + X s (SS^) 2 

+ \ Vxi S XiXi + h - c • (38) 

i 

The U(l) symmetry in the (S, y)-sector can be associated with the lepton number which 
can be spontaneously broken if S receives a VEV. The Goldstone boson related to this 
broken symmetry is known as the Majoron. Furthermore when S acquires a VEV, the 
y Xi SxiXi + h.c. term produces a Majorana mass term for x% which can be identified with the 
right handed neutrino. In this scenario the standard model neutrinos acquire a mass via a 
seesaw mechanism. The EWPT in the Majoron model has been studied in [54-56] and more 
recently in [4] including the LHC phenomenological reach. After the Higgs boson discovery 
at the LHC, the Majoron model parameters space, studied in [4], are severely constrained. 

10 EWPT and DM aspects of more general potential than Eq. (36), including softly U(l)-symmetry breaking 
terms, are studied in [50-52]. 
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4.2.1 Conformal SM + Complex Singlet Scalar + Fermions 

The EWPT phase transition in a model, closely related to the classically conformal Majoron 
model, has been recently studied in [57-59]. When requesting for 1st order phase transition, 
the model becomes rather constrained and predictive. For example, the mass of the new 
pseudo scalar ’’Majoron”, which is also a subdominant dark matter candidate, is constrained 
to be in the range ~ 1 - 2 TeV and the mass of the new CP-even scalar is ~ 100 - 200 GeV. 
The next LHC run could provide useful to investigate the model predictions [59]. 

4.3 SEyy model: SM + Complex Scalar + Fermions 

Now we will consider the EWPT in a model, which was originally build for DM phenomenol¬ 
ogy. This ”SEyy” model was introduced in [60], and features a new complex scalar field S 
charged under U(l) hypercharge, a vector like heavy electron E and a neutral fermion y, 
which is the dark matter candidate. The complex scalar has tree level interactions with 
the standard model Higgs field (Higgs portal type), whereas the dark matter particle has 
only loop-induced interactions with the standard model particles. The original motivation 
for this model was to construct an UV completion of the effective magnetic moment dark 
matter model [61] able to alleviate, at the time, the tension between the results of different 
DM direct detection experiments. The model phenomenology, including the dark matter 
and collider aspects, has been further studied in [43,60]. The possibility to secure the SM 
vacuum stability in the SEyy model context has been studied in [62] . Here we further study 
the model parameter space allowing for a strong first order EWPT. 

The difference with respect to more traditional complex Higgs portal models is that in 
our case the new state carries SM interactions 11 . 

The SEyy model Lagrangian reads [60]: 

£se XVx = £sm + X0X ~ m x xx + EilpE - m E EE - ( SE\y x + h.c.) 

+ D^D^S - m 2 s S ] S - A hs H ] HS^S - A ,(^) 2 , (39) 

where D M = <9 M — ie—Z^ + ieA^, s w and c w refer to the sine and cosine of the Weinberg angle. 
We assume that the couplings y x , A h s and A s are real. C E m stands for all the SM terms. 
The SExy x operator in Eq. 39 induces interactions between the dark matter candidate 
y and the SM fields via loop-processes and gives the name to the model. Specifically, 
quantum corrections induce an effective magnetic dipole moment A x (g 2 ) for the DM state. 
See e.g. Eq. 3 in [60] for more details. 

The scalar potential of the theory is identical to the one in Eq. (36) but S is charged 
under the SM U( 1) interactions. Depending on the values of the model parameters the model 
can feature different symmetry breaking patterns and different type EWPTs. 

If S acquires a VEV at zero temperature, the 1/(1) (hypercharge) symmetry breaks and 
the photon is massive. Thus we require the tree level potential to yield a vacuum structure 
of the type ((H), ( S )) = (v E Wi 0). Additionally we also require potential stability that leads 
to two different cases. In case (1) the S scalar does not generate a VEV at any temperature 
i.e. ( S) = 0, while in case (2) S is allowed to acquire a VEV only at some finite temperature 
i.e., ( s(T )) 7 ^ 0 at some T ^ 0 . 

n A study, considering the EWPT and collider aspects, for similar type model, including charged scalar, 
has been recently performed in [63]. 
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Figure 8: We show the results of the parameter scan for the SEyy model in case (1) producing 
a first order EWPT. The purple (blue) histograms shows the distributions of the points with 
(without) the constraints m s > 450 GeV and y x < 0.6. We have not imposed the VCs. 


4.3.1 Case (1): (s(T c )) = 0, SM like 1-loop scenario 

The tree level potential parameters satisfy the relations 

{p 2 , X,m 2 s , A s , X hs } > 0 , (40) 

so that at the tree level we have ( H) = (0, vew/V%) t and S does not generate a VEV 
at any temperature i.e. (S) = 0 and we call this case: The SM-like 1-loop scenario. The 
relevant mass terms used in the EWPT analysis are reported in App. E.l and we scanned 
the parameters within the ranges 

0.01 < Xhs < 8 , 0.01 < A s < 4 , 0.01 <y x <4. (41) 

Again we used flat prior distributions with 2000 points in the scan. In Fig. 8, we show the 
results of the scan and show that the model supports a first order phase transition occurring 
for relatively light values of the new charged scalar that can be constrained by LHC data, 
as we shall discuss below. Consistently with the Endings in the previous sections A h s tends 
to be quite large. 

4.3.2 Dark aspects and constraints 

The DM phenomenology has been studied in [43,60]. However in [43] the value y x = 2, 
motivated by the VCs, was used to constrain the model. Relaxing the VCs on y x , as we 
shall see, will improve on the model phenomenological viability. 
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We can start by noticing that the LUX experiment directly constrains the effective mag- 

2 

netic moment coupling, A x (0) ~ ■ Clearly, allowing for smaller values of y x than in [43] 

reduces also the lower limit on m s limit. Smaller values of m s are welcomed since they are 
expected to lead to a stronger first order EWPT. We also assume that mg < m x < m s which 
allows for a phenomenologically viable DM candidate. 

Following [60] we allow E to have small mixing with SM leptons (tau) such that S can 
decay into x and E, and E decays further in the SM fields. Assuming the collider limit 
for rriE > 393 GeV [60] we take m x > 400 GeV. The dominant DM annihilation channel is 
XX EE in the t — channel which occurs via the S— exchange, see Eq. 10 of Ref. [43] 


(<4vUreZ> = 


?/ 4 (m| + m 2 ) 


8n (m 


m: 


m 


2\2 ' 


mi- . , 

—y + Oly z 

A 


(42) 


For thermal (Dirac) DM the annihilation cross section is ( av re i ) ~ 4.4 x io _26 ^d for wdm > 
15 GeV [64], Using this cross section with Eq. 42 and assuming, for example rn s = 500 
GeV, m x = 420 GeV and = 400 GeV, we derive y x ~ 0.425. Such small values of y x 
and relatively light S, i.e. m s ~ 500 GeV, are now compatible with LUX that now allows 
scalar masses larger than m s ~ 90 GeV. A more refined analysis that takes into account also 
the correct thermal DM relic density leads to the constraints m s > 450 GeV and y x < 0.6. 
The purple histograms in Fig. 8 show the parameter space passing these constraints. The 
heavy electron mass does not affect the EWPT analysis unless it is heavy enough to start 
becoming Boltzmann suppressed. If this is the case, the y x dependence is suppressed in the 
S boson thermal mass and thus the transition becomes very weakly dependent on y x . 


4.3.3 Veltman conditions 

The VC for the new charged scalar reads [23] 

4A S + 2X hs + 3 g' 2 — 22/J = 0, (43) 

and for the Higgs held 

6A h + ^ 9 2 + ~^g' 2 ~ + A hs = 0 . (44) 

Using the known SM coupling values in Eq. 44 we get that A hs ~ 4.1. Inserting this in Eq. 43 
we obtain y x ~ a/2A s + 4.3. Thus, because X s > 0, for stability reasons we need to have 
y x > 2. This was the coupling used in Fig. 1 of Ref. [43] which constraints m s to be in the 
TeV range. Thus we conclude that the model supports a 1st order EWPT and a viable DM 
candidate provided the VCs are not enforced. 


4.3.4 Case (2): (s(T c )) ^ 0, tree level scenario 

Now we shall consider the second case, where the tree level potential parameters satisfy the 
conditions 

A 2 

{/i 2 , X h , —m|, A s , X hs } > 0 with 0 < A h A s < . (45) 

The latter condition is needed to ensure that the potential does not provide a VEV to S at 
zero temperature. There are also few other constraints which the parameters should satisfy 
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Figure 9: We show the finite temperature effective potential contour plot for different temper¬ 
atures: (Left) The ”EW symmetry preserving” state with a minimum at (h, s ) = (0, (s(T))) 
at T > T c ; (Center) The case with two degenerate minima (h, s ) = (h c , s c ) at T = T c \ 
(Right) The case where the EW symmetry is broken with a minimum at (h, s ) = ({h(T)), 0) 
for T < T c . This shows the occurrence of a first order EWPT in the model. 


in order for the model to support a strong 1st order EWPT. These are specified in connection 
with Eqs. 46 - 49 below. 

Once the model parameters satisfy the conditions above, it is possible that at T ^ 0, the 
system has a configuration where the Higgs has a zero VEV, but the singlet has a non-zero 
VEV i.e. the minimum of the potential is in the held space at point (h,s) = (0, (s(T))), 
where ( s(T )) ^ 0. Following [3] in this case, for the EWPT study, we can use the zero 
temperature tree level potential supplemented with the leading finite-temperature terms. 

More precisely, for a certain region of the parameter space one can directly use the tree 
level potential to establish the presence of a strongly Erst order phase transition. This 
follows from the fact that here the potential can form two degenerate minima at positions 
(h, s ) = (0, (s(T c ))) and (h, s ) = ((/i(T c )), 0) so that there is a barrier, induced by the A h. s h 2 s 2 
interaction term, between the two minima. Here the leading high temperature expansion 
terms are enough to trigger a first order phase transition between the finite temperature 
minima at (0, (s(T c ))), and ((h(T c )), 0). This potential structure for different temperatures 
is elucidated in Fig. 9. 

Intriguingly in this scenario there is a temperature range when a VEV occurs for the S but 
not for H leading to a phase in the evolution of the universe where the weak interactions are 
restored but the hyper-charge is still broken. There is, however, a new critical temperature 
above which also S vanishes and all the symmetries are restored. 

To establish the region of a strong Erst order EWPT we adopt the procedure described 
in [3] and scan the parameters using flat prior distributions in the ranges 

0.01 < X hs < 47r, 0.01 < y x < 47t, 0 < \ s , min < X s < X s ^ max < Ai r, (46) 
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Figure 10: Distributions of theSE^y model parameters in case (2) that produces a first order 
EWPT. The purple (blue) histograms show the distributions with(out) constraints m s > 450 
GeV and y x < 0.6. The VCs have not been used here to constrain the parameters. 


with the limits for A s defined as 
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s,min 


C s 

m s 


2a 2 X h 


2a\h5c s 4nA ~h$Cs^J , X SjTnax — 


X 2 

A hs 

4AV 


n s (T)/T 2 , 5c s = c s - aA s , c h = U h (T)/T 2 , 


-Vva 


s v EW 


+ | avV- c * |r ( 


(47) 

(48) 

(49) 


The parameter notation here follows closely the one used in [3] and we therefore also intruded 
the parameter a, which in our case has a value a = 1/3. The scalar thermal masses fl s (T) 
and fl h(T) are give in Eq. 95 and Eq. 96 respectively of App. E. Further, in this setup, the 
parameter ms can be written as a function of the critical temperature T c [3]. ft is possible 
to randomly choose values for T c in the range 50 GeV <T C < 250 GeV and only accept the 
values providing 1 < h c /T c < 4. Finally, A h and /i are fixed at the SM values and the scan 
is done over 5 x 10 5 points. The distributions of the model parameters producing a first 
order EWPT is given in Fig. 10. Purple (blue) histograms show the distributions with(out) 
the constraints m s > 450 GeV and y x < 0.6. These constraints were discussed in the end 
of Sec. 4.3.2. The results confirm the expectations that higher S masses tend to reduce the 
order of the EWPT if we also ask the model to provide a phenomenologically viable DM 
candidate. Performing a full one loop analysis might slightly increases the purple region. 
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4.3.5 Veltman conditions 


Enforcing the VCs of Eq. 43 and Eq. 44 requires, as in case (1), y x > 2 that is above the 
LUX bound of y x < 0.6 allowing the scalar mass to be sufficiently light for the model to 
feature a first order EWPT. 


4.4 Conformal SE%y Model 

Our final model is the classically conformal limit of the SEyy model 

U 0 = \ h {H'H) 2 + \ h9 H'H&S + \ a (S&) 2 + {SExy + h.c.) . (50) 

We assume that only the Higgs field can get a non-vanishing VEV corresponding to the 
classically conformal limit of case (1) in Sec. 4.3.1. The CW and EWPT analyses are 
similar to the one performed for the real singlet scalar Higgs portal model scenarios and the 
quantities relevant for the EWPT are shown in App. E.l. 

In Fig. 11 we show the h c /T c contour plot as function of the SEyy Yukawa coupling y x 
and the scalar S quartic self coupling A s . In the left plot, the portal coupling is fixed to 
A hs ~ 6.85 yielding the correct Higgs boson mass without VCs. The charged scalar S mass 
is therefore also predicted to have a mass m s ~ 456 GeV. If we insist on x to be a thermal 
WIMP we need to have m e < m x < m s , as argued in Sec. 4.3.2 12 . Because of the LHC 
constraints [60] we must also have 393 GeV. The DM mass is therefore constrained to 

be in a small mass window 393 GeV< m x < 456 GeV. This prediction is therefore testable 
at LHC run 2. 

4.4.1 Veltman conditions 

The Veltman conditions are again given by Eq. 43 and Eq. 44. Using the flatness condition 
A/j ~ 0 at the EW scale, and the measured values for the rest of the known SM couplings 
in the Eq. 44 we get that A hs ~ 4.84. Using this value in Eq. 43 the VC for the new scalar 
relates A s and y x . This relation is plotted as an upward turning blue curve in Fig. 11. This 
figure is our result for the ratio h c /T c as a function of y x and A s . The plot shows that the 
model can provide a strong first order phase transition when enforcing the VCs. The prize 
to pay is a too small mass for the Higgs boson which now is rrih ~ 83 GeV. 


5 Conclusions 

We analysed the EWPT of a large number of minimal extensions of the SM of particle 
interactions and their classically conformal limit. The models featured new scalars and 
fermions and we uncovered the parameter space leading to a first-order phase transition 
with(out) the Veltman-Conditions. 

When relevant we discussed their dark (matter) aspects. To accommodate simultaneously 
a first order EWPT and DM requires more involved extensions. Our results for the classically 
conformal models are summarised in Table 1. We observe that classically conformal models 
can provide at most a subdominant DM component. Allowing for non-conformal terms, 

12 We allow here hard mass terms for \ and E. These mass terms could be generated radiatively via some 
yet new unknown dynamics. 
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Figure 11: We shown the h c /T c contour plots as function of the SEyy Yukawa coupling y x 
and the scalar S quartic self coupling A s . In the left panel the portal coupling is fixed at the 
value Xhs ~ 6.85 that provides the correct Higgs boson mass without the VCs. On the right 
panel the portal coupling is fixed at the value A hs ~ 4.84, which satisfies the VC for the 
Higgs boson. The mass value of the Higgs is however not the physical one. On the upward 
turning blue curve the VC for the new scalar boson S mass is also satisfied. 


i.e. adding extra operators with dimensionful couplings, one can saturate, in certain cases, 
the full cosmic DM density. In Table 2 we also summarise the results for few non-conformal 
model examples studied here. In the main text the reader will find further discussions of the 
results for the models summarised in the tables as well as for other minimal extensions not 
reported in the tables. 

A positive feature is that all these extensions must feature new particles lighter than a 
TeV making these models testable at the upcoming LHC run two experiments. 
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Conformal model 

Higgs mass 

1 st order EWPT 

DM 

SM 

- 

+ 

- 

SM + VC 

- 

+ 

- 

RS; (s) = 0 

+ 

+ 

subdom. S [34] 

RS; (s) ^ 0 * 

+ 

+ 

- 

RS + VC; (s) ^ 0 * 

- 

+ 

- 

RS + Wf; (s) = 0 

+ 

+ 

- 

RS + Wf + VC; (s) = 0 

+ 

+ 

- 

RS + Wf; (s) ± 0 

+ 

+ 

subdom. Mf 

RS + Wf + VC; (s) ^ 0 

+ 

+ 

- 

CS; (s) = 0 

+ 

+ 

subdom. S [34] 

CS + Df; (s) = 0 (SEyy) 

+ 

+ 

- 

CS + Df + VC; (s) = 0 (SEyy) 

- 

+ 

- 

CS + Mf; (s) = 0 [59] 

+ 

+ 

subdom. PS [59] 


Tabic 1: Summary of the different conformal models that can provide the observed Higgs 
boson mass, 1st order EWPT and DM. The Higgs portal models can feature Real Singlet 
(RS) with(out) a Weyl fcrmion(s) (Wf), or a Complex Scalar (CS) with(out) either a Dirac 
fermion(s) (Df) or Majorana fcrmion(s) (Mf). We will also consider the cases with and 
without adding the Veltman Conditions (VC). In these models the DM particle is identified 
as a scalar (S), a pseudoscalar (PS) or a Majorana fermion. The model, RS + VC; (s) = 0, 
is not included in the above list because the associated one-loop Higgs mass squared is 
negative. In the last column we also refer to the works where the results for the DM part 
were obtained. The results for the model in the last row are taken from [59]. Finally the 
symbol # means that the model suffers from the cosmic domain wall problem. 


Non-conformal Model 

Higgs mass 

1 st order EWPT 

DM 

SM 

+ 

- 

- 

SM + VC * 

+ 

- 

- 

RS + Wf + VC; (s) = 0 (PNC) 

+ 

+ 

+ Mf [43] 

CS + Df; (s) = 0 (SEyy) 

+ 

+ 

+ Df 

CS + Df + VC; (s) = 0 (SEyy) 

- 

+ 

+ Df 


Table 2: Summary of the different non-conformal models that we have considered in this 
work and that can yield the observed Higgs boson mass, a 1st order EWPT and possibly a 
DM candidate providing the full cosmic DM density. We follow the notation introduced in 
Table 1, with the difference that here # means that the model does not lead to the observed 
top quark mass. 
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Appendices A-F 


We summarise in the following appendices the effective potential and other relevant 
quantities needed to determine the order of the EWPT for the various models investigated 
in the main text. In the last appendix the relevant equations needed to determine the DM 
relic density for the PNC model are summarised. 

A Effective Potential 

The full one-loop ring improved finite temperature effective potential, used in our perturba¬ 
tive study of the order of the EWPT, can be schematically written in the form 

V ES (h, s, T) = V 0 (h, s) + V 1 {h, s) + V T (h, s, T ) + 5V(h, s) . (51) 

Here V 0 denotes the tree level potential, whereas V\ and Vt are the one-loop zero and finite 
temperature parts of the potential respectively, and 5V is the counter term. We give details 
of these different terms below. We will write the effective potential as a function of the 
background fields ( h,s ), that are related to the Higgs doublet H via H = (0, h/\/2) T 1 and 
to the new scalar field S via S = s/a , where the normalization factor a = 1 for real scalar 
and a = \/2 for the complex scalar. 

We specify the tree level potential Vo for each model separately at the beginning of each 
section in the main text. To be precise, we implicitly take that V 0 (h, s ) = Vsm(^) + V 0 (h, s ), 
where Vsm is the SM part of the potential, except for the pure scalar part of the potential, 
which is given in Vo. Also the non-SM contributions to the potential are defined in Vo- 
The one-loop zero temperature contribution to the effective potential is the Coleman- 
Weinberg one 


U(M = X> 


rii 


mf(h,s) 
647r 2 


log 


' rn%{h,s) 

Q 2 



(52) 


where the sum runs over the bosons (+) and fermions (—) and n* counts the internal degrees 
of freedom of each species i. The background field dependent tree level masses for each field 
are denoted by rrii(h,s), Q is the renormalization scale, C % is a constant equal to 5/6 for 
gauge bosons and 3/2 for scalar bosons and fermions. The relevant degrees of freedom and 
specific mass terms rrii(h, s ) are given for each model separately in App. C, D and E. 

The counter term part of the potential can be written in the form 


SV(h, s ) = 5 -^h 2 + °-^h A + + + t 


SX, 


5 X h 


5 m 2 


SX* 


(53) 


where the counter-term parameters 5/d 2 ,5Xh,5Xh s ,S7nl,5X s , and 5Vq are fixed by the renor¬ 
malization conditions. It is more practical to define these conditions for the classically 
conformal and non-conformal models separately. 

Among the non-conformal SM extensions we have studied in detail the ”SEyy” model, see 
Sec. 4.3. For case (1) of this model presented in Sec. 4.3.1, the renormalization conditions 
are similar to that of the SM i.e. only 5/i 2 and SXh are relevant and need to be fixed by 
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requiring that the minimum of the effective potential and the Higgs boson mass value are 
held at their tree level values, that is 


d_ 

dh 


(V 1 (h,s) + 5V(h,s)) 


h=VEW 

s=0 


0 and 


d 2 
< 911 2 


(V 1 (h,s) + <5V(h,s)) 


h=VEW 

s=0 


0 . 


In case (2), studied in Sec. 4.3.4, because we are using only the tree level potential in 
the EWPT study (with leading finite T contributions), we do not need to determine the 
counter-terms. The renormalization conditions for the conformal models are specified in 
Sec. A.l. 

The one-loop finite temperature part reads 

rjn 4 POO 

Vr{h,s,T) = J2 ±rii ^ J d ^ lo S 

+ [™>K h > s ) - m i(h, s, T)] , (55) 


1 =F exp — \jy 2 + rrif(h, s)/T 2 


(54) 


where again (+) refers to bosons and (—) to fermions. The last term including the Kronecker 
delta function takes into account the usual ring improvement (daisy diagram resummation) 
for bosonic (B) degrees of freedom, where m i(h, s, T ) stands for the thermally corrected boson 
mass. To speed up the calculations we used an approximation for the integral in Eq. (54) [42] . 
This approximation consists in smoothly joining the high and low temperature expansion of 
the integral and is summarized in App. B. 


A.l Classically conformal case 

In the classically conformal cases the EW symmetry breaking needs to be triggered by loop 
effects i.e. via dimensional transmutation a la Coleman-Weinberg (CW) [9]. For this to 
happen the tree level scalar potential needs to have a flat direction. The flat direction 
conditions, for generic scalar potentials has been provided in [65]. Particularly, we follow 
the flat direction and CW analysis done in [23] for the models studied here, which relies 
on [9,65]. We shall briefly summarise the core points here and end with the renormalization 
conditions for the classically conformal case. 

Let us start with classically conformal tree level scalar potential 

V = + c.t. , (56) 

where the \jki are renormalized couplings and c.t. refers to counter terms. In principle the 
model can (and will) feature also fermionic and vectorial contributions but for our purposes 
to find the flat directions in the pure scalar potential it is enough to consider the previous 
form of the potential and take that \jki ~ g 2 ~ y 2 1, where g and y refers to possible 
gauge and Yukawa couplings included in the model. For the CW mechanism to be operative 
the model scalar quartic couplings need to be the such that there is a flat direction at tree 
level. Particularly the couplings \ijki should be renormalized such that at the renormalization 
scale Q the potential has indeed a flat direction. To this end, let rq be a unit vector in the 
scalar field space, then to find the flat direction we need to solve 

min (Ajjfc/'Uj'Uj'Ufc'U;) mui=i 0. (^0 
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If a solution Ui = n t can be found, then 0* = defines a flat direction in tree level potential. 
Thus the couplings should be renormalized such that 

A ijki(Q)rnn j n k ni « 0 . (58) 

Now we shall write the potential using a classical background field such that (j) c i = n^i. 
For our analysis it is convenient to use polar coordinates in the scalar held space. The tree 
level potential acquires the form 

V(h,s) o = V{r,u) 0 = r 4 f(uj), (59) 

where the function f(oj) can be read from the original potential V(h,s )o once substituting 
h = rcoscu and s = r shun. As usual, the minimum can be found by studying the first 
and second derivatives of the potential. Specifically the model parameters satisfying the 
condition f'{uj) = 0 define the direction (uj) of the non-trivial minimum as discussed in 
Sec. 3.2.3. Now the unit vector defining the hat direction is n = (cos (a?), sin (a;)). And the 
hat direction problem in our case takes the form 


min ( XijkiUiUjUkUi ) \ UiUi =i = mmf(uj)\ UiUi=l = 0 . (60) 

As we are using the polar coordinates the background held is now (j) c i = r, and the one loop 
potential V\ (h, s ) can be written in the form 


Vi (<f>d,u) = 
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Q 2 
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= A4^ + B(f>% log 


Q 2 
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and B ^ £ ± n, W?(w) 


6Air 2 


(61) 

(62) 


(63) 


and we have decomposed all the mass terms in the form mf((j) c i,uj) = Wf {uj)^ using the 
fact that in the classically conformal case all the mass terms are directly proportional to (f) c i. 
The quantity W?(uj) thus depends only on the model couplings since in the end also uj will 
depend only on the scalar quartic couplings. 

The full zero temperature potential now reads V((J) c i,uj ) = V 0 ((p c i,uj ) + Vi (jp c i, uj) + 
5V((pci,w)- And finally, the renormalization conditions for quartic couplings, fixing the 
counter terms, is defined via relation 

= A ijki{Q)ninjn k ni = 0 , (64) 

m? (<j>d ,td)=Q 2 e( _25 / 6+c i) 



with the renormalization scale chosen to be Q ~ vew■ The factor e~ 25 / 6 is extracted for 
convenience and follows the renormalization choice made in [65]. For completeness, let us 
also give the relation for the scalar (Higgs) mass in this case. The non-trivial extremum of 
the potential Eq. 62 is at 


log 



1 A 

2 ~ B ’ 


(65) 
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and if this corresponds to a ground state, the scalar fluctuation along the flat direction has 
a positive mass squared 


m cw — , 

which defines the CW scalar boson mass. 


( 66 ) 


B Thermal integral 

Here we give the approximation we used for the thermal integral, Eq. 54 that smoothly joins 
the high and low temperature expansions of the thermal integral. This approximation has 
been used previously in [42]. The thermal integral in Eq. 54 is of the form 


IbA 0 ) = ± / dr/r/ 2 log 
■Jo 

where the upper sings refer to bosons and the lower to fermions and a = m 2 (/r, s)/T 2 . The 
high temperature expansions ( rrii/T <C 1) of the integral, Eq. 67, for the bosons and fermions 
can be written as [66] 

4 2 2 

t ii. . IT 7T 7T 3 a ,, , x N , _ _ 

J f0) = Q a2 - 32 0 [o M - c B ), (68) 

r rr. . 7vr 4 7T 2 a 2 , , . . 

Ip (a) = + —a + — (log(a) ~ c F ), (69) 

where c E — 3/2 — 2^ E + 2 log(47r) and c E — 3/2 — 2j e + 2 log(7r) and 'y E ~ 0.5772. The low 
temperature expansion (m*/T 1) of the integral is the same both for bosons and fermions 

and reads [48] 


1 Te -V^+~ a 


(67) 


= _ yf a3/4e • ( 70 ) 

The approximated expression Ip F for the thermal integral, I e ,f reads 

^B, F ( a ) = tB,F(a)lB,F( a ) + (1 — JB,FIB,F ~ IbA 0 ) i (71) 

where the fitting functions for bosons and fermions, that join the high and low temperature 
expressions, are £g(a) = e~( a / 6 - 3 "> and ig(a) = e _< - a//3 ' 25 ^ respectively [42], The relative 
error in this approximation at worst is less than 5% in the range 0 < a < 20. Finally, the 
approximated finite temperature part of the effective potential is 


Vf(h,s,T) 



' B F 

'22n i l£(m?(h, s)/T 2 ) + ^n^(m 2 (h, s)/T 2 ) 


Yl n ‘i2i [ m <' s ) _ m ?( ft . s . r )] • 


(72) 

(73) 


where n l counts the internal degrees of freedom and in the second line we have added the 
thermal ring improvement for the bosonic degrees of freedom. 
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C Mass Terms for the SM 


Here we give the relevant mass terms used in the one-loop effective potential, Eq. 51, when 
studying the EWPT in the SM. These terms are also relevant for the other models considered 
in this work. The background held dependent masses, discussed in connection with Eq. 51, 
are of the form mj(/i, s) = rrii(h) and rrii(h, s,T) = mi(h,T), because the only background 
held in this case is the Higgs mean held h. Now, the well known h dependent tree level mass 
terms in the SM are 


m 2 ho (h ) = 3A h h 2 - /i 2 , m 2 G (h ) = X h h 2 - /r 2 , (74) 

m w>T (h) = m W:L (h) = | h, m z ,T{h) = m z>L (h ) = + 9 h, (75) 

m t (h) = ^J=h, (76) 

where m,h 0 and me are the Higgs boson and the Goldstone boson masses respectively. We 
have divided the gauge boson masses mw,i and rriz.i into transverse (i = T) and longitudinal 
{i = L) components for later use. Finally, m t is the top quark mass. As mentioned before, we 
only take into account the top quark contribution from the SM fermions. For completeness, 
we write down also the Goldstone boson terms, even though we do not consider them in 
the analysis, as discussed in the end of Sec. 1. The internal degrees of freedom for different 
species are 


nh 0 = 1, no = 3 (three Goldstones), 

n w,T = 4 (polarization x particle/antip.), nw,L = 2, n z ,T = 2, n Zt L = 1, 
n t = 12 (spin x particle/antip. x color), 


(77) 

(78) 

(79) 


where we have used the same subscript notation as for the masses above. The scalar boson 
thermal masses are 


■<y.T) = + nf"(T), (80) 

m %(h,T) = m%(h) + nf'tT), (81) 


where the thermal self energy (at momentum p = 0 limit) for the Goldstone and for Higgs 
boson is the same 

n : s h M (T) = ^ [3 g 2 + g' 2 + 8 A h + Ay 2 ] T 2 . (82) 

Finally, the longitudinal gauge boson thermal masses follows by diagonalizing the mass 
matrix 


M 2 L (fc,T) 


/ 

9 2 

0 

0 

0 

\ 


( 

s 2 

0 

0 

0 

\ 


0 

g 2 

0 

0 


h 2 

T + 


0 

g 2 

0 

0 



0 

0 

9 2 

-99' 



0 

0 

s 2 

0 


V 

0 

0 

-99' 

g' 2 

) 


V 

0 

0 

0 

g' 2 

) 



(83) 


where the two upper rows refers to W-boson contributions, and the other two to the Z-boson- 
photon mixing contributions. The gauge boson transverse components do not get thermal 
corrections at the one-loop level, and thus the transverse thermal masses are equal to the 
tree level masses i.e. m 2 T (/r, T) = m 2 T (h ) for i = W, Z, 7 with m 2 T (h ):s given in Eq. (75) 
except for photon, which is = 0 . 
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D 


Mass Terms in the Extension with A Real Singlet 
Scalar 


D.l Case (1) and non-conformal scenario 

In this case the singlet scalar does not acquire a non-zero VEV and the EWPT analysis is 
similar to the SM case. The difference with respect to the SM case resides in the corrections 
to the SM scalar boson mass terms and the addition of the singlet scalar effective potential. 
The latter is effectively obtained, i.e. the part contributing to the EWPT, by adding the 
singlet in the sums in Eq. 52 and Eq. 54 with a tree level mass term m 2 Q {h) = m 2 s + A hs^ 2 , 
and one internal degree of freedom, i.e. n s = 1. The scalar boson thermal masses are now 


(h,T) 

= ml o (h)+Il h (T), 

(84) 

(h, T ) 

= m 2 Jh) + MT), 

(85) 

(h, T ) 

= m 2 G (h ) + n h (T), 

(86) 


where the thermal self energies (at momentum p = 0 limit) are 


I MT) 

n,CO 


n s h M {r) + —t 2 , 


A hs 

IT 


4 


A hs, 

12 

T 2 


(87) 

( 88 ) 


with n^ M (T) given in Eq. 82. Notice also that there is no mixing between h and s in this 
case. Thus, in the end, the physical Higgs held is directly related to h. 

In the classically conformal case, i.e. case (1), the parameters are as above except that 
p 2 and m| are set to zero. 


D.2 Case (2) 

Here the singlet scalar acquires a non-zero VEV, the mass terms given above are therefore 
partially modified. The tree level scalar mass terms are obtained by diagonalizing the mass 
matrix 


M ho,s 0 


{h, s ) = 


A/> 'Mi T A/j S s 2X hs hs 


2A hs hs X s 3s 2 + X hs h 2 

and the scalar thermal masses by diagonalizing the matrix 


Mft 0)S0 (h, s, T) — M^ Q So (h, s) + 


MT) 

0 


0 

MT) 


(89) 


(90) 


where the thermal self energies are given in Eq. 87 and Eq. 88. All the other mass terms are 
like in the SM, with p 2 and m| set to zero. 

As discussed in the main text and in App. A.l, for the CW and EWPT analysis it 
is convenient to use polar coordinates (r,uj) in the scalar held space. We get the mass 
parameters in these coordinates by substituting h —y rcoscu, and, s —> rsincu. 


35 




D.3 Mass terms in PNC Model 

D.3.1 Case (1) 

In this case the singlet scalar does not acquire a VEV and therefore the relevant terms are 
like in case (1) considered in Sec. D.l, except for the extra contribution following from the 
singlet Weyl fermion. Specifically, the only modification is in the singlet scalar thermal self 
energy (at momentum p = 0 limit), which now reads 


n,(T) 


hs , A s 1 9 

T + 4 +3^ 


rp2 


(91) 


D.3.2 Case (2) 

Here also the singlet scalar receives a VEV, in this case the mass parameters are like in case 
(2) considered in Sec. D.2 above, except for the (Weyl) Majorana fermion x that receives a 
tree level mass m x (s ) = 2 y x s, and should be taken into account in effective potential, with 
n x = 2 . Of course we need to take into account also of the singlet scalar thermal self energy 
IL a (T) given in Eq. 91. 


E Mass Terms in SEyy Model 


E.l Case (1) and conformal scenario 

In this case the scalar S does not acquire a VEV at any temperature. Thus the EWPT follows 
closely the SM and the real singlet Higgs portal scenarios, in cases (1), considered above. 
As in the case in Sec. D.l, here the only background held is the h held, and the background 
held dependent masses are of the form rrii(h,s) = rrii(h) and rrii(h, s,T) = mi(h,T). All 
terms related to the SM particles are just like in the Sec. C above, except for the mass terms 
which we shall refine here. Thus just like in the case of SM + real singlet, addition to the 
SM contributions to the 1-loop effective potential, there is now extra piece coming from the 
S scalar, with tree level mass term m 2 (h ) = m 2 s + with internal degrees of freedom 

n s = 2 for complex scalar. The scalar boson thermal masses are 


m l(h,T) 

= m 2 h (h) + U h (T ), 

(92) 

m l(h,T) 

= m 2 (h) + n s (T), 

(93) 

m Uh,T) 

= m 2 G (h)+ U h (T), 

(94) 


where the thermal self energies (at momentum p = 0 limit) are 


n*(T) 

n.CO 


nr CO + t '\ 

V + h + I,« +1,/ 

6 3 4 y 12 y 


T 2 


(95) 

(96) 


with U s h M (T) like in the SM in Eq. 82. Notice that there is no mixing between h and s in this 
case. Thus the physical Higgs field is associated directly with h. Finally, the longitudinal 
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gauge boson thermal masses we get by diagonalizing the mass matrix 


M l(h,T) = 


( 

9 2 

0 

0 

0 

\ 


( 

9 2 

0 

0 

0 


0 

9 2 

0 

0 


h 2 


0 

9 2 

0 

0 


0 

0 

9 2 

- 99 ' 


T + 


0 

0 

9 2 

0 

V 

0 

0 

- 99 ' 

g ’ 2 

J 


V 

0 

0 

0 

g ,2 { 1 + 


\ 


11 2 
6 ’ 


(97) 


where the two upper rows refers to W-boson contributions, where as the two lower to the 
Z-boson-photon mixing contributions. The new contributions to the B-field thermal mass 
follows from the new charged fields: from the scalar S and from the vector like electron E. 
The transverse thermal masses are as in the SM case i.e. mf T (h, T) = m 2 T {h) for i = W, Z, 7 
with m 2 T {h )\s given in Eq. (75) except for photon which is m 7 ,r(/r) = 0. 

I 11 the classically conformal scenario, the mass terms are as above, expect that the di¬ 
mension full scalar Lagrange mass parameters, /j 2 and m|, are set to zero everywhere. 


E.2 Case (2) 

Following the approach of [3], for the tree level strong 1st order EWPT calculations, in 
addition to the tree level potential, we need only the one-loop scalar thermal masses. These 
are given already in Eq. 95 and Eq. 96 above. Further, the tree level model parameters 
should satisfy the relations given in Eq. 45 related to Eqs. 46 - 49 in Sec. 4.3.4. 

For completeness, we also give the mass terms in a form that can be used for a full one 
loop EWPT calculation. 

Because we have two background fields, h and s, the mass terms depend on both, 
i.e. rrii(h,s) and rrii(h,s,T). We have also mixing between the h and s fields here and 
thus we write the mass terms in matrix form. We deduce the tree level scalar masses by 
diagonalizing the scalar mass matrix 


/\/f 2 fi \ — I ^h,3h 2 — /i 2 H—f^s 2 

h ' sS ^ ' A Hs hs A s 3s 2 - ml + ^h 2 


4 . 2 


The Goldstone boson masses are 

A h h 2 - jj 2 + ^s 2 


M^(h,s) = 


0 


0 


A s s 2 


m 


l + %Or 2 


(98) 


(99) 


where the first row, associated with G, stands for all the three SM Goldstone boson masses, 
whereas the lower row, refers to a, that is the mass of the Goldstone boson associated with 
the broken U(l) symmetry. This occurs because the electrically charged scalar S acquires a 
background field dependent VEV. The Thermal masses of the scalars are 


M l f [h,s,T) = Ml,{h,s) + 


n,,iT:i 0 \ 

0 n .iT ) ) ■ 


( 100 ) 


where Ll/ l (T) and L1 S (T) are as in Eq. 95 and Eq. 96 respectively. The tree level gauge 
boson masses are as in the case (1) above, except that the B-field has now a mass term for 
non-zero s values. Obviously, in this case, the photon has a tree-level .s-induced background 
mass term. The Z-boson-photon mass mixing term, which holds for both longitudinal and 
transverse components, at tree level is now 


M 2 


z, 7 


(M) = 


9 2 ~9Sf \ 4 

- 99 ' g ' 2 J 4 + 


0 

„/2 
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0 

0 


( 101 ) 






And finally the longitudinal gauge boson thermal masses are 


M 2 L (h,s,T) = M 2 L {h,T) 


/ 0 0 0 0 \ 

0 0 0 0 

0 0 0 0 

V o 0 0 g' 2 ) 


( 102 ) 


where M \{h,T) is as for case (1) of Eq. (97). 


F Relic density and annihilation cross sections in PNC 
Model 


Here we give the relevant equations needed to evaluate the y relic density for the PNC model, 
case(2), studied at the end of Sec. 3.4.4. 

To obtain an estimate for the y relic density we solve the usual DM Boltzmann (Lee- 
Weinberg) equation using the classic approximate method. Here, the freeze-out temperature 
Tf is solved iteratively via the equation 


x f 


m y , m Y 

rr log 2A 


/ 45 77T 2 , 

-— ((TUrel), 

2 g* x f 


(103) 


where m v \ is the Planck mass and g*{T) is the number of effective energy degrees of freedom 
for which we use the standard model values. For the thermally averaged annihilation cross 
section we use the expression [68] 

(<TtW> = Iftn^K/r) Li V s ~ 4m » gl ds, (104) 

where K\ and K 2 are the modihed Bessel functions of the second kind and the annihilation 
cross section to different final states i, a = cr xx ^, are given below in Eqs. 106,108,114 (and 
below 115) and 116. The relic density for the y particle is obtained from equation 


Qh 2 


1.07 x 10 9 x f 
\fg*m p i(av ie i) ‘ 


(105) 


The above approximate solution for the relic density does not work well near resonances, 
but for our purposes it yields a reasonable overall estimate of the DM density. We use it 
for the full range assumed by the y masses that can yield a first order EWPT. A more 
accurate solution for the relic density can be obtain by solving the Lee-Weinberg equation 
numerically, however we do not expect this to affect our conclusions, that is, y can be at 
most a subdominant DM component. 

Now we summarise the y annihilation cross sections into different final states including 
SM fermions, Higgs and gauge bosons and the <f> scalar. 

The annihilation cross section into fermionic final states reads 


a 


xx - H / 


y 2 sin 2 uj cos 2 to m 2 Xf 


^ nV EW 


+ 


s mj, + s — m 2 ^ + iT 


sM /, ( 106 ) 
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where the factor Xf — 1 for leptons and for quarks 

X f =N c f l + ^(3l og (!ti)+. 


37T 


(107) 


which takes into account the relevant one loop QCD correction for light quarks in the range 
s < 90 2 [67]. The color factor N c = 3 and we use the value a s = 0.12 for the strong coupling. 
The velocity factor is /% = \Jl — 4 rrif/s. 

The annihilation cross section into gauge bosons is 






— Ib 


?/ 2 sin 2 u> cos 2 unri% 

27Tv‘^ 
l EW 


+ 


s — + iT $171$ s — + iY 


x 


PxPb j—j- - —y + 3 


4 m A B 


m 


B 


(108) 


(109) 


where B refers to W or Z bosons and the factor fw = 1 and fz = 1/2 respectively. 

The trilinear scalar couplings, needed when calculating the tree level annihilation cross 
sections into scalar final states, are defined as 


A= vXh cos 3 uj + uXhs sin u cos 2 uj + vXh s sin 2 c o cos u ; + uX s sin 3 uj 
A$$$ = — vXh sin 3 u + uXhs sin 2 uj cos uj — vXh s sin uj cos 2 uj + uX s cos 3 uj 
X $00 = 3nA s sin 2 uj cos uj — 3 vXh sin uj cos 2 u + 

Xh s {u cos uj — 2u sin uj cos uj — v sin uj + 2v sin uj cos uj) 

A^$$ = 3uX s sin uj cos 2 uj + 3vXh sin 2 uj cos uj + 

A hs [v cos uj — 2u sin u j cos uj + u sin uj — v sin u j cos uj) , 

where v = uew and u = wew tana;. 

The annihilation cross section into the scalar (j) (Higgs) final states is 


( 110 ) 

( 111 ) 

( 112 ) 

(113) 


a xx-^H ~ 


y; 


647 ts 2 /3? 


2 yl sin 4 u 


2 Logf (m/ + 4m 2 (s - 4m 2 )) 


sfi<j>Px I ~b 


s — 2 m | 

y\ sin 4 uj Logf (—2ni0 + 8m 2 + s) + 

stt^x (% m i - 4 ( 8m x + s ) m l + 64m 4 + s 2 - s 2 /3|/3 2 ) 


2yt sin 4 uj 


s2 PlP 2 x ~ ( s - 2m l) ‘ 


+ 


Logf (-2m| + 8m 2 + s) J + 

16 m x y x sin 2 u j (Logf (2m 2 — 8m 2 + s) + 2s /3<t,f3 x ) x 
3 sin uj (m 2 — s) X 


(m 2 -s) + 


16s (4m 2 - s) 


cos uj (■m 


s ) A 


<I>( t>(j) 



(ml - s) + m|T| 
+ 


9 sin 2 ujXI h 


cos 2 ujX% h 


s) + m^T 2 (m% - s) + m% T 2 


+ 


6 cos u sin uX^^X^ ((m 2 - s) (m% - s) + m^m^Y^Y^) 


m: 


s ) + m l T l) U m l - s ) +m%Y% 


(114) 
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where 


Logf = log 


$Wx + !) - 2m l 


s {faPx ~ !) + 2m 0 


(115) 


The annihilation cross section into heavier $ scalars, <r xx _>.$$, is obtained by interchanging 
silica cosca and 0 •<-)•$ in the above cross section, cr xx _^, given in Eq. 114. 

Finally, the annihilation cross section into (j) + $ final states reads 


1 




647 ts 2 /3?, 


-4jd sin 2 (ca) cos 2 (ca) 


(m 2 + m 2 - s) {-P x y/\~f + ml + m 2 - s) (/3 x ^X~ f + m% + m 2 - s) 


x 


LogF (-s (m| - 16m 2 + 3m^) + 2 (m| - 4m 2 ) (2m| + 4m 2 + m 2 ) + s 2 ) x 

(^ A / - i m l + m l - s ) 2 ) - 2 PxV\f ( m l + m l ~ s ) x 

- 3) m^ + /3 2 (m| - s ) 2 + ml (-2/1 2 (m| + s) - 2m| + 16m 2 + 3s) + 


where 


2 4 on 4 2 ' 

sm $ — m $ — 32m x — s 


16m x 7/ x sin(ca) cos(ca) (^2/3 x -\/Xf — LogF (— 8 m 2 + 2m| + s) j x 

sin(ca) (s - mj) cos(ca) (m| - s) 

T 2 m 2 + (m 2 - s ) 2 T|m| + (m 2 - s)" 


+ 


^ 4 . 4 , cos 2 (ca) 


+ 


A l<«, sin 2 (ca) 


8/3y\ZXfyl (s — 4m 2 ) . 

\^r|m| + (m| — s) r 2 m 2 + (m 2 -s) 

2A$^A^$$ sin(ca) cos(ca) (r$r^m$m^ + (m| - s) (m| - s)) 


+ 


r|m| + (m| - s)‘ 


ml - s 


(116) 


LogF 

A / 


= log 


m| + m 2 - s 


- s -/^x\/ A / 
m| - s) 2 - 2m 2 (m| + s) + m^. 


(117) 

(118) 
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